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^ ABSTRACT 

We present a catalogue of structural parameters for 8814 galaxies in the 25 fields of the 
HST/ACS Coma Treasury Survey. Parameters from Sersic fits to the two-dimensional sur- 
face brightness distributions are given for all galaxies from our published Coma photomet- 
ric catalogue with mean effective surface brightness brighter than 26.0 mag arcsec"^ and 
brighter than 24.5 mag (equivalent to absolute magnitude - 10.5), as given by the fits, all in 
F814W(AB). The sample comprises a mixture of Coma members and background objects; 
424 galaxies have redshifts and of these 163 are confirmed members. The fits were carried 
• • out using both the Gim2d and Galfit codes. We provide the following parameters: Galaxy 

^ ^ ID, RA, DEC, the total corrected automatic magnitude from the photometric catalogue, the 

total magnitude of the model (F814Wab), the geometric mean effective radius Re, the mean 
surface brightness within the effective radius (/i)e, the Sersic index n, the ellipticity and the 
source position angle. The selection limits of the catalogue and the errors listed for the Sersic 
parameters come from extensive simulations of the fitting process using synthetic galaxy mod- 
els. The agreement between Gim2d and Galfit parameters is sensitive to details of the fit- 
ting procedure; for the settings employed here the agreement is excellent over the range of 
parameters covered in the catalogue. We define and present two goodness-of-fit indices which 
quantify the degree to which the image can be approximated by a Sersic model with concen- 
tric, coaxial elliptical isophotes; such indices may be used to objectively select galaxies with 
more complex structures such as bulge-disk, bars or nuclear components. 

We make the catalog available in electronic format at Astro- WISE and MAST. 

Key words: galaxies: clusters: individual: Coma; galaxies: elliptical and lenticular, cD; galax- 
ies: dwarf; galaxies: fundamental parameters; galaxies: evolution 



Based on observations made with the NASA/ESA Hubble Space Tele- 
scope, obtained at the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc. under 
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1 INTRODUCTION 

Surface brightness distributions are a vital tool in our understand- 
ing of galaxies. Since the pioneering work of Reynolds (1913) 
and Hubble (1930) on elliptical galaxies it has become common 
to fit the radial surface brightness distributions to functions hav- 
ing a small number of parameters, which include a scale length, a 
characteristic surface brightness, and one or two further parameters 
which describe the structure of the surface brightness profile. The 
most commonly used fitting function is that of Sersic (1963, 1968), 
whose function includes as special cases both the i?^/^ law of de 
Vaucouleurs (1948) and the exponential surface brightness distri- 
bution which is characteristic of disk galaxies (Patterson 1940; de 
Vaucouleurs 1957, 1959; Freeman 1970): 

I{R) = le exp{-6 [(i?/i?e)'/" - l] }, (1) 

where I{R) is the specific intensity at distance R from the centre. 
Re is the radius enclosing half the galaxy light, le is the specific in- 
tensity at i^e, ^ is the Sersic index or concentration index (Trujillo 
et al. 2001), and h ^ 1.9992 x n - 0.3271 (Capaccioh 1989). 

The Sersic function provides a good model for ellipticals, gi- 
ants showing values of n ^ 4, intermediate luminosity ellipti- 
cals n ^ 2 — 4 and dwarfs n ^ 1 — 2 (Caon, Capaccioli & 
D'Onofrio 1993; Graham et al. 1996; Graham & Guzman 2003). 
Bulges of disk galaxies are also well fit by the Sersic model (An- 
dredakis, Peletier & Balcells 1995) with indices n ^ 0.5 — 4 (Bal- 
cells, Graham & Peletier 2007; Graham & Worley 2008). For disk 
galaxies, pure Sersic fits often yield poor approximations to the 
entire galaxy surface brightness distribution, due to the presence 
of bulges, bars, spirals, outer disk truncations (e.g. van der Kruit 
& Searle 1981a,b) and anti-truncations (Erwin et al. 2005). How- 
ever, classifying galaxies into early (n > 2.5) and late (n < 2.5) 
types on the basis of single Sersic fits to the entire galaxy has be- 
come standard practice (e.g. van der Wei et al. 2008), especially 
in samples with limited image depth and spatial resolution which 
prevent more complex modeling. This practice fails when the sam- 
ples include lower luminosity dwarf elliptical galaxies. The relia- 
bility of such fits may be calibrated by performing single Sersic 
fits to nearby, well resolved galaxies. Needed for the interpretation 
of such single Sersic fits is a parameter that quantifies the degree 
to which the true surface brightness distribution deviates from the 
Sersic model. 

The HST/ACS Treasury Survey of the Coma cluster was pre- 
sented in Carter et al. (2008, Paper I). Although the survey was 
originally planned to cover 740 arcmin^ of the Coma cluster field, 
the final areal coverage is 274 arcmin^ in the F475W and F814W 
bands, mostly in the core region, owing to the ACS failure in 2007 
January. Still, with the exquisite quality and depth of the imaging 
and the large number of spectroscopic redshifts known for galaxies 
in this field (Colless & Dunn 1996; Mobasher et al. 2001; Marzke 
et al. 2010; Chiboucas et al. 2010) this survey allows studies of the 
structure of large samples of cluster members to an unprecedented 
depth. The photometric catalogue from the HST/ACS images was 
presented in Hammer et al. (2010, Paper II; see Sect. [2]). 

This paper presents a structural analysis of sources selected 
from a structural analysis of the sources from the Paper II photo- 
metric catalogue, based on two-dimensional single Sersic fits. Of 
the ~75,000 objects in that catalog, we provide Sersic parameters 



NASA contract NAS 5-26555. These observations are associated with pro- 
gramme GO 10861. 



for 8814 galaxies that are located both inside the cluster and in the 
background; the selection function is explained in Sect. |6.1| We 
present standard Sersic parameters as well as two goodness-of-fit 
indices, that provide a quantitative measure of the degree to which 
the galaxy surface brightness distribution deviates from a Sersic 
model with concentric, co-axial elliptical isophotes (Sect. [5]). These 
indices can be used to identify those galaxy images which allow for 
additional components, such as outer disks, nuclear components or 
bars. Given the complexity of the structural analysis process, we 
focus this paper on the presentation of the analysis techniques and 
of the catalogue. We defer scientific analysis to future papers. The 
structural parameters presented here can be used: 

• To study the cosmological evolution of galaxy sizes and 
shapes by using the Coma cluster as a local reference sample. 

• To quantify the faint end of global scaling relations, such 
as size-(surface brightness) diagrams and the Fundamental Plane, 
revealing how dwarf elliptical galaxies do or do not unite with 
brighter ellipticals. 

• To study the correlation between the structural parameters and 
the photometric masses of elliptical and lenticular galaxies, which 
could be used in cluster membership studies (Trentham et al. 2010). 

Our results in Coma can be compared with the lower density 
Virgo and Fornax cluster environments where targeted HST/ACS 
surveys provide structural information at higher physical resolu- 
tion for smaller samples of galaxies (Ferrarese et al. 2006; Cote et 
al. 2007). The Coma data set may be also used in conjunction with 
HST surveys at higher redshift to study the evolution of the struc- 
tural properties of galaxies. STAGES (Gray et al. 2009) is a survey 
of the supercluster Abell 901/2 at a redshift of 0.165. Amongst an 
extensive multi-wavelength dataset, ACS images have been used 
for Sersic fits to a large sample of galaxies in the STAGES region. 
GEMS (Rix et al. 2004) is an ACS survey of a 900 arcmin^ region 
within the Extended Chandra Deep Field South region. Although it 
is a field rather than cluster survey, it provides a useful evolution 
benchmark at redshifts approaching z = 1. HST has been used to 
study the structural properties of galaxies in higher redshift clus- 
ters, where there is a suggestion of size evolution by up to a factor 
four (e.g. Trujillo et al. 2006; Strazzullo et al. 2010). 

There are currently a number of codes capable of performing 
two-dimensional Sersic model fits to the surface brightness distri- 
bution of galaxies. Two extensively used are Gim2d (Simard et 
al. 2002) and Galfit (Peng et al. 2002). Both codes work by min- 
imising a merit function, and produce similar outputs, but their in- 
ner workings differ in a number of ways, such as the minimisation 
technique: Gim2d uses the Metropolis algorithm, whereas Gal- 
fit uses the Levenberg-Marquardt algorithm. GALFIT offers prac- 
tical advantages, such as higher execution speed and the ability to 
simultaneously fit several targets. But because each code has its 
own merits, we carried out the fits using both codes, and present 
both results. In order not to bias the comparison, two teams worked 
largely independently, one with Gim2d, and another with Galfit. 
While some details differ, e.g. in the parameter ranges explored in 
the Montecarlo simulations, there was enough coordination to en- 
sure the results would be comparable to each other. We show in 
Sect.|7]that the agreement is very good. 

Gim2d and Galfit were compared by Haussler et al. (2007), 
who concluded that both could produce similar results, but warned 
against a systematic underestimate of both the total luminosity and 
effective radius of the Gim2d output for lower surface brightness 
sources. We were able to reproduce their findings; in Sect. |3] we 
present an approach which successfully overcomes these biases. 
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2]describes the input data 
and §[4] describe how the 



The paper is structured as follows, 
from which our catalogue is derived. §[3]; 

Gim2d and Galfit analysis runs on the Coma data were setup. §[5] 
introduces two additional parameters calculated from the residuals 
of the data from the models which describe how well the Sersic 
model fits the data. §[6]presents the final structural parameters cat- 
alogue and the criteria for inclusion in the catalogue. §[?] presents 
a comparison between the results obtained with two codes on the 
Coma data. In §|8]we present our conclusions and describe the next 
steps in the analysis of the galaxies and their surface brightness 
distributions. In Appendix[A|we compare our results with those of 
Haussler et al (2007). 

Throughout the paper we assume the distance to Coma of 100 
Mpc, corresponding to a distance modulus ofm — M = 35.0 (see 
Paper I). All magnitudes are in the AB system. In this paper, with 
a few exceptions, we express the surface brightness in terms of 
the mean effective surface brightness (/i)e, i.e., the mean surface 
brightness enclosed within Rq. Graham & Driver (2005) show that 
the relationship between (/i)e and the effective surface brightness 
//e (surface brightness at Re) is given by: 



(M)e = Me - 2.51og[/(n)], 
where: 



h poo np>^ 

, e-^x'"-Mx = ^r(2n) 



(2) 



(3) 



and r(2n) is the complete gamma function. Total magnitude m is 
related to (/i)e by the simple relation: 

m= (/i)e-2.51og(27ri?^), (4) 



2 DATA 

2.1 HST/ACS images 

The survey design and reduction of the images are described in de- 
tail in Paper I, so only a summary will be provided here. The obser- 
vations (program GO 10861) were taken between 2006 November 
and 2007 January with the HST/ACS camera (2x4096x2048 pix- 
els, 0.05 arcsec/pixel, 3 arcsec interchip gap). A total of 25 visits 
were completed before the failure of ACS. Most fields (19/25) are 
located within 0.5 Mpc of the cluster centre (the full list of survey 
fields is given in Table 2 of Paper I). The remaining six fields are in 
the South- West extension of Coma. Two HST orbits were devoted 
to each pointing. A four-position dither pattern was used for each 
of the F475W and F814W images, with total integration times of 
2560 s and 1400 s, respectively. The dither pattern allowed us to 
fill the ACS inter-chip gap, albeit with lower S/N. Total exposure 
times were lower for some visits due to dither positions that failed 
to acquire guide stars. Final exposure times are given in Table 5 of 
Paper I. 

Data reduction was carried out with a dedicated pipeline. It in- 
cluded the combination of individual images with the Multi-Drizzle 
software (Koekemoer et al. 2003), which yields combined images 
resampled onto a rectified (but original sky orientation) output 
frame with 0.0495 arcsec/pixel. Cosmic rays were removed during 
the multi-drizzle process and also using LACOSMIC (van Dokkum 
2001). These processed images, together with an initial source cat- 
alogue comprised the first data release (DRl), 2008 Augus{^ 

^ MAST (archive.stsci.edu/prepds/coma/) and Astro-WISE (www.astro- 
wise.org/projects/COMALS/) 



2.2 Catalogues of the Coma Data Release 2 

The second data release (DR2), available at the same web sites 
as DRl, includes improvements in alignment between F814W and 
F475W images, better astrometry, aperture corrections to the SEx- 
TRACTOR photometry, and photometry of sources that project onto 
large galaxies. Details of the data processing and description of the 
DR2 photometric catalogues, including the S EXTRACTOR (version 
2.5; Bertin & Arnouts 1996) configuration parameters employed in 
the catalogue generation, are given in Paper II. The DR2 catalogues 
contain ^73,000 sources. Based on Monte-Carlo simulations, the 
80% completeness limit for point sources in the DR2 catalogues is 
27.8 mag in F475W and 26.8 mag in F814W. 

The DR2 images and catalogues are the basis for the structural 
analysis done with Galfit, whereas Gim2d fitting was performed 
on the DRl images and catalogue. This difference represents no 
problems. Comparison of the SEXTRACTOR catalogues from the 
two releases shows that 99% of the detected sources match, with 
75% of the additional catalogue objects in DRl being in Visit 03, 
which lacked two of the four dither positions. These additional ob- 
jects would in any case be fainter and smaller than the catalogue 
limits of the current paper (Sect. |6.1| l. We refer to the DR2 cata- 
logue as the "photometric catalogue" throughout this paper. 



2.3 PSF 

The Point Spread Function (PSF) is a key ingredient of any signal- 
to-noise weighted analysis of the morphological and structural 
properties of galaxies. The PSF of the Wide Field Channel (WFC) 
of the ACS has been extensively studied. Jee et al. (2007) and 
Rhodes et al. (2007) offer different suites for creating ACS PSFs for 
a variety of observing conditions. The HST Intrument Science Re- 
port 0306 (Krist 2003) presents a detailed study of the variation of 
the PSF across the WFC chips. The PSF of the WFC depends both 
on time and position on the chip. The TinyTim program (Krist 
1993) takes advantage of this empirical knowledge, and creates ar- 
tificial PSFs for a large variety of observing conditions and HST 
instruments. 

We created a grid of ACS PSFs using TinyTim. These were 
then combined using the code DrizzlyTim, by Luc Simard, 
kindly made available to us by the author. 

DrizzlyTim calculates the location of the original PSFs in 
the calibrated flat fielded individual exposures, using the same mul- 
tidrizzling parameters and shift file as used to produce the science 
images. DrizzlyTim then invokes TinyTim to create the required 
PSFs in the calibrated, flat fielded set of coordinates. The PSFs are 
created with an oversampling of 5, and assuming a 6500 K black 
body as a representation of the object spectrum. This is appropriate 
for the E and SO galaxies in our sample. DrizzlyTim then places 
these PSFs into blank frames, with the same size and header pa- 
rameters as those of the real flat fielded individual exposures. These 
frames are then coadded, using again the same multidrizzle param- 
eters as those used to manufacture the final science images. The 
final step is to apply the Charge Diffusion Kernel to these newly 
created PSFs in the calibrated, geometrical distortion corrected im- 
ages. We followed this process to create a grid of PSFs with the 
sampling of the original images. One PSF was created every 150 
pixels in x and y directions, and the PSF imagelets created were 31 
pixels on a side. The typical FWHM of the PSFs created was 2.0 
pixels, with at most a 20% variation across the field. As the time 
period over which the images were obtained was short, we did not 
allow for temporal variation of the PSF. When fitting real galaxies. 
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Figure 1. Four examples of the PSFs created by Drizzly Tim. These PSFs 
were used by Gim2d and Galfit to fit the real sources in the science im- 
ages. The image in each panel is 50 pixels (2.0 arcseconds) square. 



both Gim2d and Galfit were instructed to use the nearest PSF to 
that object. 

Fig.[T]shows a set of DrizzlyTim PSFs, with their respective 
FWHM. 



3 GIM2D FITTING 

This section describes the Gim2d setup used to fit the surface 
brightness distributions of the detected galaxies. As will also be 
done in Sect. |4] for Galfit, we build upon the experience of the 
GEMS collaboration who analyse galaxies from HST/ACS imag- 
ing, and perform extensive simulations to compare systematics in 
Gim2d and Galfit (Haussler et al. 2007). From our own analysis 
of the Haussler et al. paper, shown in Appendix [A] we conclude 
that mask creation is critical for Gim2d, and we present a new 
prescription for creating the masks that Gim2d requires for object 
detection and fitting. In addition, this section presents MonteCarlo 
simulations used to assign statistical errors to the Gim2d best-fit 
parameters. 



3.1 Object masks for Gim2d 

Our own simulations, and those of Haussler et al. (2007), show that 
Gim2d can miss a substantial fraction of light from faint sources 
(a comparison of our simulations to those of Haussler et al. is given 
in Appendix [A|. Our analysis concludes that this problem origi- 
nates in the use of S EXTRACTOR segmentation masks as input for 
Gim2d. When Gim2d is instructed to infer its initial guesses from 
S Extractor magnitude and size parameters (D0INIT=YES), 
and is allowed to refine the sky level estimate obtained from S EX- 
TRACTOR (DOBKG=YES), the use of SEXTRACTOR segmentation 
masks leads to systematically fainter and smaller solutions. Setting 



DOBKG and DO IN IT to NO (i.e. sky background fixed from SEX- 
TRACTOR, and intitial estimates of other parameters taken from the 
parameter file gal.mdpar) fixes the systematic error in total magni- 
tude and Re, but at the expense of an increase scatter in the so- 
lutions, and a dramatic increase in the convergence time. When 
DOBKG and DOINIT are both set to YES, Gim2d is left free to 
automatically decide which section of the global parameter space 
to explore. Gim2d does this using the mask image it has been pro- 
vided. It first estimates the sky using the pixels designated as sky 
pixels in the mask image that lie a fixed number (which we set to 
10) of pixels away from the target mask. It then subtracts this es- 
timate of the sky value from the input image, and derives initial 
estimates of the total flux, inclination angle, ellipticity and effec- 
tive radius. In this important step, Gim2d calculates the total flux 
and the effective radius of the target using pixels which, according 
to the mask it has been given, belong to the ISOAREA of that ob- 
ject. Gim2d will then explore the range in parameter space from 
counts to twice the sky subtracted flux within the mask, and from 
pixels to twice the effective radius of the set of pixels within the 
mask. This factor of two is hard coded into Gim2d and can not be 
tweaked. 

Therefore if Gim2d is set with both DOBKG=YES, 
DOINIT=YES and fed with a segmentation image from SEX- 
TRACTOR, it will only explore magnitudes between (MAG_ISO - 
0.75) and infinity. However, the SEXTRACTOR simulations pre- 
sented in Paper II show that for the fainter sources the real mag- 
nitudes can be off from MAGJ^UTO by up to two magnitudes. A 
very similar statement could be made for the effective radius. In 
this case, the range in linear size explored is from pix to twice the 
Gim2d initial estimate, which is built from the two dimensional 
Kron radius. This effective radius estimate was found to be differ- 
ent from the true effective radius by up to a factor of 5 for the less 
luminous sources. 

This clearly indicates that the masks Gim2d is provided have 
to be enlarged, if they are to be a faithful representation of the real 
extent of the targets. Instead of using the standard segmentation im- 
age created by SEXTRACTOR, we build a customized mask for each 
object, using the information from the SEXTRACTOR catalogue for 
the whole frame and the knowledge of the noise properties of the 
images. This mask image is constructed separately for each par- 
ticular object, since Gim2d treats target and background sources 
differently. Any given object is represented by one aperture when 
acting as the target, and is represented by another smaller aperture 
when being considered a background source possibly affecting the 
fit of a different object. The properties of the proposed masks, to- 
gether with the practical steps required to create them, are summa- 
rized in Appendix [B] 

Fig. [2] shows four examples of galaxies, together with their 
associated masks. Black pixels belong to the target object, white 
pixels belong to other sources, and grey pixels are sky. The first ob- 
ject is a low surface brightness source that was however detected by 
SEXTRACTOR. The second and fourth objects are spiral galaxies. 
In all cases, the FOVs are given in the image insets. The Kron-like 
apertures which were adopted for the targets, presented as black 
pixels on the mask images, are clearly more extended than the 
visible flux from the target galaxy, and are typically much larger 
than the isophotal apertures used for neighboring objects which are 
shown as white pixels. The apertures for the background sources 
can overlap the aperture of the target object, which is not possible 
when using the segmentation images produced by SEXTRACTOR, 
when the background objects can potentially interfere with the abil- 
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Figure 2. Examples of the masks used by Gim2d. The "Kron-like" aper- 
tures devised for the 4 target objects are much larger than the aper- 
tures of their neighboring sources, which have an area equal to their 

ISOAREA_IMAGE. 



ity of Gim2d to measure regions outside the ISOAREA_IMAGE of 
the target. 

3.2 Noise model 

A simple noise model is used for the fits. When Gim2d is not given 
a specific noise image, it builds an internal weight map based upon 
the rms of the background (cTbkg) and the effective gain, which 
is a function of the effective exposure time of the science frames. 
This noise model is a transaction of the usual CCD uncertainty 
equation. In our case, abkg is taken directly from the S EXTRAC- 
TOR catalogue, although it is later recalculated by Gim2d in most 
cases. The effective exposure time is read from the header of the 
HST/ACS frames. 



With DOBKG set to YES, Gim2d refines the sky value given 
by S Extractor, and obtains a better estimate of abkg- This is 
then used to construct the internal weight map. The sky pixels in- 
volved in this calculation were at least 10 pixels away from those 
pixels determined to belong to the target object. The refined sky 
value is the median value of at least 30 pixels, applying a 5cr clip- 
ping thresholding scheme. The rms of the background obtained by 
Gim2d in this way is generally in excellent agreement with the 
cTbkg initially estimated by SEXTRACTOR, both in the mean and 
variance. In cases where insufficient numbers of pixels were avail- 
able to estimate the background rms inside the imagelet (owing 
to crowded fields), GIM2D defaulted to the user-supplied rms that 
was estimated by SExtractor. 

3.3 Final Gim2d configuration 

DO IN IT was set to NO. The lower limit of the total flux was set 
to 0, while the upper limit was set to 10 times the automatic aper- 
ture flux. The Bulge- to-Total fraction was set to 1. The effective ra- 
dius ranged from 0.0 to 10.0 times the effective radius estimate ob- 
tained for a pure Sersic model of n = 2.25 of the same MUOBS and 
FLUX_RADIUS. The ellipticity and position angle were allowed to 
search their whole ranges. The X and Y drifts were permitted to 
range from —10 to 10 pixels, and the residual sky value was al- 
lowed to go from —0.25 to 0.25. The Sersic index n was allowed 
to range from 0.25 to 10.0. Thus, each object had an individualized 
Gim2d configuration file leading to all objects being fit by a single 
Sersic profile. The Metropolis temperatures are adequate for the ex- 
plored ranges, and expected typical changes in each iteration. In all 
cases, saturated pixels were rejected from the fits. The Metropolis 
algorithm is given 400 iterations to cool off after achieving conver- 
gence (see Simard et al. 2002 for more technical detail). 

3.4 Errors on the parameters and depth of the survey 

Although Gim2d produces, together with its results, a set of confi- 
dence intervals for the fitted parameters, the error estimates repre- 
sent only the scale upon which the Figure-of-Merit that Gim2d 
uses is expected to vary. Therefore these confidence intervals 
merely reflect how constrained the fit is. A more realistic and mean- 
ingful error analysis needs to investigate the extent to which the 
minimum of the Figure-of-Merit can drift in its parameter space. 

A modest number of Monte Carlo Gim2d simulations was 
run. The purpose of these simulations is twofold. The first is to 
be able to ascribe realistic statistical errors to the fits produced by 
Gim2d. The second purpose is to assess the surface brightness limit 
beyond which it will not be possible to recover reliable structural 
parameters. 

10,000 model images were created using GALIMAGE within 
IRAF.FUZZY. In this step, a Moffat (1969) PSF representative of 
the average properties of several non-saturated stars was used to 
degrade the galaxy models, mknoise was then used to add appro- 
priate Poissonian noise to this model. This blurred image was then 
added to a real ACS image which therefore provides the readout 
noise, and the bulk of the error correlation that is typical of ACS 
drizzled images. The selected canvas images correspond to visits 1, 
15, 78, and 90 (see Paper I). SExtractor and Gim2d were run, 
with the same parameters, weight images and flag images as those 
used to create the SExtractor catalogue and the same experi- 
mental setup described above, using the Moffat PSF as the convo- 
lution kernel. 



6 Carlos Hoyos et al. 



The mean effective surface brightnesses of the models ran 
from 19.0 to 27.0, with effective radii distributed randomly in 
log Re between 2.0 and 60.0 pixels. Sersic indices were randomly 
distributed between 0.5 and 4.5. Ellipticities were randomly dis- 
tributed from 0.0 to 0.8 and position angles were unconstrained. Of 
the total number of model galaxies created, 10,000 were both de- 
tected by S Extractor and successfully analyzed by Gim2d; the 
analysis presented in this subsection deals with these 10,000 fake 
sources. The remaining sources were either not detected by S Ex- 
tractor or fell in problematic areas of the image such as the CCD 
edges and were thus rejected. 

The first step in the analysis of these simulations is presented 
m Fig. [3] This fi gure shows the magnitude residuals against the ef- 
fective radius residuals, and the effective radius residuals against 
the Sersic index residuals. These relations are presented in differ- 
ent panels, according to the input n as indicated in the figure. The 
points with an input (/i)e brighter than 22.0 are highlighted in blue. 

As expected, the residuals in total magnitude anticorrelate 
with those in effective radius: where Gim2d yields a higher than 
expected luminosity, it also yields a higher than expected effective 
radius. This occurs for all Sersic indices. The oblique, green line 
shows the error correlation that would preserve the mean effective 
surface brightness within one effective radius. The figure shows 
that Gim2d introduces a surface brightness bias: overluminous so- 
lutions have fainter (/i)e. The slope of this covariance is found to 
be Sersic index-dependent. The contraction or expansion of the fit- 
ted functions with respect to the input parameters is not arbitrary, 
it depends on the real profile of the underlying source being fitted. 
Although these observations merely confirm an expected behaviour 
given the experiment being run, they are important because they al- 
low us to use these simulations to assess the errors on the output 
parameters. 

Comparable behaviour can be seen in the right panels of 
Fig. [3] If Gim2d finds Re higher than the real value, the output n 
is also higher than the real value. This happens because n is essen- 
tially determined by the pixels in the object wings, which are much 
more abundant and have higher weighting due to the lower Poisson 
noise, thus Gim2d has to increase the model power in these wings. 

In addition, the fact that the Sersic index has such a large im- 
pact on the behaviour of the residuals implies that this parameter 
has to be included in our recipe for error assessment (see also Mar- 
leau & Simard 1998, their figure 1 1). 

Fig. |4] shows the magnitude, effective radus and Sersic in- 
dex residuals against both input and output mean effective surface 
brightness for the 10000 simulations. This figure is divided into 
four panels: for Sersic indices of 0.5 < n < 1.5, 1.5 < n < 2.5, 
2.5 < n < 3.5, and 3.5 < n < 8.0, respectively. The black 
lines with error bars give the run of the median value of the resid- 
uals and a robust estimate of the vertical la scatter of the residuals 
around this latter median value. These median values of the resid- 
uals and their associated la scatter have been calculated in equally 
spaced bins in input (left panels) and output (right panels) mean ef- 
fective surface brightness. For clarity, the vertical la scatter values 
are plotted separately in Fig. [5] 

From the data presented in Figs.|4]and[5]we find that, for (/i)e 
^24.0, the median values of A log(n), A \og{Re) and A mag (the 
systematic errors) are always less than 0.05, 0.04, and 0.08, respec- 
tively. These values are only weakly dependent upon n, with the 
highest Sersic index bin having smaller systematic errors (less than 
0.05, 0.03, and 0.05 respectively). For brighter surface brightness 
((/i)e ^ 21.0) these differences are lower than 0.02, 0.03, and 0.04. 
The widths of the distributions, which we associate with the non- 
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-4.43 


0.14 
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-4.93 


0.16 


logn(G2D/Model) 


Sersic index 




a 




0.5 < n < 1.5 


-3.81 


0.12 




1.5 < n< 2.5 


-2.81 


0.08 




2.5 < n< 3.5 


-2.74 


0.07 




3.5 < n< 8.0 


-3.11 


0.09 



Table 1. Table of coefficients required to use Equation |5] to estimate the 
statistical errors on the total magnitude, Re and n, for various ranges of 
output n, for the Gim2d fits. 



systematic error in the recovery of the input values, increase toward 
lower surface brightness, and are, at the faintest limit, 0.12 dex, 
0.15 dex, and 0.25 mag, for log(n), \og{Re) and total magnitude, 
respectively. These values and Figures indicate that, with the use of 
the tailored masks, Gim2d is indeed able to recover the input pa- 
rameters accurately. The use of these customized masks and the use 
of the individualized search in the parameter space allows Gim2d 
to have a better understanding of the galaxy flux and size. This nat- 
urally leads to a better fit, free from the systematic errors that were 
detected by Haussler et al. (2007) and confirmed in Appendix [A|. 

Figs. |4] and [5] indicate that it is reasonable to use the output 
Sersic index and (/i)e to derive realistic error estimates to the total 
magnitudes, effective radii and Sersic indices. Given the modest 
number of simulations, we adopt a simple two-parameter approach 
based upon the output (/i)e and n. A single straight line of the form: 

loga = a X (/i)e,out + /3, (5) 

is fit to the robust la vertical scatter around the median shown in 
the right panels of each quadrant of Fig. [5] Although this functional 
form is expected for the magnitude uncertainties, it is also used 
for the uncertainties in Re and Sersic index for simplicity. Table [T] 
shows the best-fit coefficients for these fits. 

To assign meaningful and realistic statistical errors to any 
Gim2d measurement, we first calculate the output (/i)e. Next we 
evaluate the linear functional form given above using the coeffi- 
cients found in Table [T| The final uncertainty in the parameter of 
interest, as given in the structural catalogue, is the antilogarithm of 
the result. 

Finally, the run of the errors with (/i)e and n, given in Fig.|5] 
and derived from the coefficients in Table[T] allow us to infer a lim- 
iting output (/i)e beyond which Gim2d will not be able to success- 
fully recover the true parameters. We choose an operational limit of 
an uncertainty of 0.25 mag. This limit corresponds to S/N < 5.0, 
and given the coefficients in the table, the corresponding limiting 
(/i)e is 24.5 mag arcsec"^. Since the magnitude is the first moment 
of the light distribution of any object, it will not be possible to re- 
liably recover the remaining structural parameters, which would be 
higher moments, from lower surface brightness objects. 
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1814W,AB(G2D"Mdl) log 10(Re(G2D/Mdl)) 



Figure 3. Effective radius residuals vs. magnitude residuals (left panels) and Sersic index residuals vs. effective radius residuals (right panels) as a function of 
the input Sersic index and {/i}e. The residuals are defined as the ratio of the output value from Gim2d to the input model value. From top to bottom, the panels 
include sources with input Sersic index 0.5 < n < 1.5, 1.5 < n < 2.5, 2.5 < n < 3.5, and 3.5 < n < 8.0, respectively. The blue dots represent points 
with input {fi)e brighter than 22.0, and the black dots represent the whole pool of models created. The green line shows the location of the points should the 
fitting process preserve the mean effective surface brightnesses of the models. Deviations from this line indicate that Gim2d was not able to accurately retrieve 
the value of (/i}e- Also, the slope of the clouds is correlated with the intrinsic profile of the model being fit. For galaxies with large n, it is somewhat more 
difficult to reproduce the parameters of the input model. The right set of panels shows how the fits to the models expand or compress, depending on whether 
the fit overestimates or underestimates Re . 



4 GALFIT FITTING 

This section describes the Galfit (version 2.0.3c) setup used to 
fit the program galaxies as well as the simulations carried out for 
error assessment. Nearly all galaxies included in the SEXTRAC- 
TOR catalogue presented in Paper II were fit, except for the sources 
that were originally buried in the extended haloes of large galaxies. 



Galfit is capable of fitting multiple galaxies simultaneously. Be- 
cause its minimisation algorithm is based on a gradient method 
it is significantly faster than Gim2d. However the algorithm is sus- 
ceptible to getting stuck in a local minimum. For the fit to converge 
quickly to the correct values it is essential that the initial values of 
the parameters are as close to the real solution as possible. 

Most initial parameters that we use for fitting are based on the 
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Figure 4. Magnitude, effective radius and Sersic index residuals for the Gim2d simulations. This figure is divided into four quadrants, each with six panels. 
Top left quadrant (a) shows the residuals for simulations of galaxies with input 0.5 < n < 1.5; top right quadrant (b) simulations with input 1.5 < n < 2.5, 
bottom left quadrant (c) simulations with input 2.5 < n < 3.5, and bottom right quadrant (d) simulations with input 3.5 < n < 8. In each quadrant, the top 
two panels show the residuals in Sersic index as a function of input (/i}e (left) and output (/i}e (right) panels. The middle two panels show the residuals in 
Re and the bottom two the residuals in magnitude, again against input and output {/i}e. The lines with vertical error bars show the run of the median value of 
the residuals; the error bars are 1.5 times the inter-quartile width of the vertical distribution. The colour coding shows the two dimensional histogram of the 
density of the underlying points, normalised along the vertical axis only. The lowest level (black) has a density >1% of the maximum, and the highest level 
(red) is >50% of the maximum. Intermediate shades are at 5%, 10%, 15%, 20%, 30% and 40%. 



S Extractor catalogue. When fitting large numbers of galaxies 
any manual intervention is extremely time consuming. Therefore, 
we decided to make use of ASTRO-WIS^ which provides a facility 
for the structural analysis of large datasets. 



www.astro-wise.org 



4.1 Galfit setup in Astro-wise 

Astro -WISE is an information system and environment for large 
imaging datasets, up to the Petabyte regime, with multiple users 
around Europe. In ASTRO-WISE one can archive raw data, calibrate 
data and perform scientific analysis storing all results. Valentijn et 
al. (2007) provide a technical description of the information sys- 
tem, and Sikkema (2009) describes the data reduction pipeline. 

The scientific analysis components in ASTRO-WISE include. 
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Figure 5. Amplitude of the error bars for the Gim2d simulations plotted in Fig.|4] The figure layout is identical to that of Fig.|4] In each panel, the right 
panels provide a direct measure of the (/i}e dependence of the non-systematic error in each of the model parameters A log(n), A \og(Re) and Amag. 



among others, routines for source extraction (using SEXTRAC- 
TOR), variability analysis, photometric redshifts and galaxy sur- 
face photometry fitting using Galfit. The Galfit implemen- 
tation within Astro -WISE automatically produces the specified 
postage stamps of the sources, runs Galfit itself and stores the 
configuration and results for all sources in a database. ASTRO- 
WISE enables full backward-chaining of data lineage in general. 
This means that model image, residual image and customized in- 
spection plots can be created anytime upon request. 

We use the Galfit setup of Haussler et al. (2007) as a starting 
point for our fitting scheme. Here we present a detailed description 
of the fitting process. 

First, Astro-wise makes a postage stamp of each source. 
The size chosen is slightly larger than that of Haussler et al., and 



is determined from the SEXTRACTOR image size measurements, 
such that 



size = 4 • A_IMAGE • KRON_RADIUS. 



(6) 



Next, a sigma image is created from the Inverse Variance Map. This 
sigma image is modified to take into account Poisson noise from the 
sources as well. 

Nearby sources are masked according to the SEXTRACTOR 
segmentation image. In line with the GEMS results and encouraged 
by the results with Gim2d we expand the masks for nearby sources 
by using elliptical apertures with semi-major axis 4xA_IMAGE 
and with ellipticity and position angle as determined by SEXTRAC- 
TOR. Sources for which this mask overlaps with the mask of the 
main source are fitted by Galfit together with the main source. 
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The implementation of the Sersic profile in Galfit has 8 free 
parameters. Of these, we leave the diskiness/boxiness parameter 
fixed so that all isophotes describe perfect ellipses. All other pa- 
rameters are left free. In addition to this, we leave the sky free, 
although we do not allow for any gradient in the sky. We do, how- 
ever, constrain n to the interval [0.5,8.0]. Gradient based fitting 
methods do require an initial guess for all parameters. Except for 
the Sersic index, which we initialize as n = 1.5 for all sources, 
initial guesses for each source are based on parameters from the 
SEXTRACTOR catalogue: for Re we use FLUX_RADIUS [3] , for 
total magnitude we use MAG_I SO. The axis ratio and position angle 
are initialized from ELLIPTICITY and THETA_IMAGE. 

Astro-wise uses these parameters to write a configuration 
file for Galfit. In the case of mulitobject fitting, we determine the 
input parameters for the secondary objects in the same way, with 
the exception that we keep the position of the source fixed if its 
centre is outside the postage stamp. 

4.2 Shot-noise Simulations 

Similarly to what was done for Gim2d (Sect. |3.4| ), an extended 
set of simulations were performed. The simulations serve three 
main purposes. First, they allow us to test our Galfit setup, by 
identifying biases in the fits. Second, they allow us to infer realis- 
tic errors of the output parameters. Like Gim2d and other fitting 
codes, Galfit tends to underestimate errors on the fitted parame- 
ters (cf. Haussler et al. 2007). Our simulations allow us to assign er- 
rors to fitted parameters which are more realistic than the standard 
Galfit errors. (Our errors are still lower limits because images of 
real galaxies deviate from the perfect Sersic model with concen- 
tric, coaxial isophotes.) Finally, the simulations allow us to define 
limits for the minimum signal-to-noise required for reasonable fits. 
The simulations were not designed to test the performance of Gal- 
fit in crowded areas, which has already been extensively discussed 
by Haussler et al. (2007). 

Galfit is wrapped in Astro-wise using the python lan- 
guage which allows for straightforward customization and script 
writing for a specific science case. We adapted the python code in 
Astro -WISE to create simulated galaxies, insert them into images, 
create source lists and then to run Galfit on them. As is the case 
with Gim2d, Galfit's ability to correctly fit a given galaxy varies 
with the intrinsic parameters of the fitted galaxy. Hence, to assign 
errors to the fit parameters of real galaxies we require the results of 
a large number of simulated galaxies with similar output parame- 
ters. 

In our approach we created a mock catalogue of 200,000 
galaxies. The parameter ranges used are listed in Table [2] Each pa- 
rameter samples the given range, either uniformly or uniformly in 
the log as indicated in the Table. The parameter ranges were cho- 
sen so that the distributions of output parameters bracket the dis- 
tributions found in the data. When generating these parameters, we 
avoided the edges of the frame and applied a hard cut-off in fie to 
avoid any detection problems with SEXTRACTOR. 

After each model galaxy had been fitted by Galfit, the distri- 
butions of the differences output minus input were binned in order 
to determine the variation of the errors with key output parameters 
such as magnitude, surface brightness and Sersic index. We chose 
six bins in magnitude, ten bins in log(i?e), two bins in ellipticity 
and five bins in n (600 bins in total). To minimise the uncertainty 
on the errors, one would like to have as many galaxies per bin as 
possible. Our 200,000 models yield ~330 per bin, so that, even 
though for some galaxies the output bin will be different from the 



parameter 


range 


#bins 


X (pix) 


400.0... 4000.0 




y (pix) 


400.0... 4000.0 




mag 


20.0... 25.0 


6 


Re 


2.0. . . 60.0 (log) 


10 


n 


0.5... 6.0 (log) 


5 


Me 


<25.5 




ell 


0...0.8 


2 


pos 


0.0... 180.0 





Table 2. Input simulation parameters of Galfit single Sersic galaxies. All 
quantities are cast uniformly in the given range, except Re and n, which 
are cast uniformly in logarithmic space. The number of bins denotes into 
how many bins the range was divided for the final error assesment. In the 
Galfit simulations the effective surface brightness rather than the mean 
effective surface brightness was used to define the bins. 

input bin and certain output bins will be more sparsely populated 
than others, the relative uncertainty on the errors (assuming they 
are Gaussian) is always less than 10%. 

Galfit itself was used to generate the artificial galaxy mod- 
els. Although one might argue that using GALFIT to make the two- 
dimensional images to which itself it should fit models is doubt- 
ful, we stress that Galfit has been tested extensively and that, in 
our opinion, it is doing at least a better job than IRAF.ARTDATA, 
which does not sufficiently oversample in the centres of galaxies. 
Our simulation setup takes into account convolution of the model 
galaxies with a DrizzlyTim PSF. Before injecting them into real 
ACS observations, Poisson noise was added to these galaxies. To 
avoid any crowding, we used only 100 models per ACS-frame, so 
that we ended up with 2000 frames, each with 100 artificial galax- 
ies on top of the ^^2500 sources already present. Simulated galaxies 
were injected into Visit 90, because this frame is relatively empty 
and on does not suffer from any missing dithers. 

On each frame with simulated galaxies we ran SEXTRACTOR 
using the same configuration as was used for the real data (see Pa- 
per II). We associated our list of simulated sources with the sources 
detected by SEXTRACTOR by demanding that they be at most 14 
pixels away from the closest source in the catalogue. A small frac- 
tion (^1%) of sources were not detected by SEXTRACTOR. In a 
small number of cases, SEXTRACTOR can be confused by proxim- 
ity to or even blending with a source already present in the frame. 

Results of the simulations are shown in Figs. [6] and [T] As in 
Fig.|4] Fig. [6] shows residuals of the logarithm of the Sersic index 
log(n), the logarithm of the effective radius log(i?e), and the to- 
tal magnitude, against input (left) and output (right) mean effective 
surface brightness {p)e. In general, the distributions of the output 
parameters around the mean have non-Gaussian, extended wings. 
Hence, a standard RMS error does not allow for a straightfor- 
ward interpretation. Rather than determining the RMS of a outlier- 
clipped sample, we use a 95% confidence interval determined from 
the interquartile range per surface brightness or magnitude bin. The 
symmetrized intervals are used as errorbars in the plots in Fig. [6] 
and plotted again in Fig.[7] 

The results are excellent, and, overall, similar to those ob- 
tained with Gim2d (Sect [34]). Up to {p)e = 24.0, the median 
differences, or systematic errors, are below 0.04 dex, 0.02 dex, and 
0.04 mag, in log(n), \og{Re) and total magnitude, respectively, 
except for the highest Sersic index bin where the differences at 
{p)e — 24, are always lower than 0.06 dex, 0.06 dex or 0.07 mag. 
The slight bias pattern that appears for {p)e > 24.0 and high Sersic 
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Figure 6. Magnitude, effective radius and Sersic index residuals for the Galfit simulations. This figure is organized as Fig.|4] Magnitude, effective radius 
and Sersic index residuals for the Gim2d simulations. This figure is divided into four quadrants, each with six panels. Top left quadrant (a) shows the residuals 
for simulations of galaxies with input 0.5 < n < 1.5; top right quadrant (b) simulations with input 1.5 < n < 2.5, bottom left quadrant (c) simulations with 
input 2.5 < n < 3.5, and bottom right quadrant (d) simulations with input 3.5 < n < 6. In each quadrant, the top two panels show the residuals in Sersic 
index as a function of input (/i}e (left) and output {/i}e (right) panels). The middle two panels show the residuals in Re and the bottom two the residuals in 
magnitude, again against input and output (/i)e. The lines with vertical error bars show the run of the median value of the residuals. As in Figure gthe error 
bars are given by 1.5 times the interquartile range. The colour coding shows the two dimensional histogram of the density of the underlying points, normalised 
along the vertical axis only. The lowest level (black) has a density >1% of the maximum, and the highest level (red) is >50% of the maximum. Intermediate 
shades are at 5%, 10%, 15%, 20%, 30% and 40%. 



indices, (Fig. [6j:,d), is an boundary effect of the simulation setup. 
Output models are brighter, and have larger Re, than the input val- 
ues due to the fact that input models only reach (/i)e < 24.0. The 
region with output {/j)e > 24.0 is only populated with models for 
which Galfit has found a solution with fainter (/i)e- Because the 
errors in Re and (/x)e are coupled, these models must have positive 
Re residuals, as observed in the middle-right panels of Fig.|6j:,d). 

The run of non- systematic errors with input and output mean 
effective surface brightness (Fig. [7]) shows similar behaviour to 



those of the Gim2d errors (Fig. [5]). For a given surface brightness, 
Galfit errors tend to be slightly smaller than Gim2d errors but 
the differences are not meaningful, given that Gim2d values are 
more uncertain owing to the lower number of Gim2d simulations. 

In Table [3] we present the parameters needed to estimate the 
uncertainties of the GALFIT output parameters, as was previously 
done in Table [T]for GIM2D. This gives the coefficients for the fits 
of Equation[5] The best-fit relations were then applied to extract the 
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Figure 7. Amplitude of the error bars for the Galfit simulations plotted in Fig.|6] The figure layout is identical to that of Fig.[6| In each panel, the right 
panels provide a direct measure of the (/i}e dependence of the non-systematic error in each of the model parameters A log(n), A log(i?e) and Amag. 



error estimates for each of the Sersic parameters that are given in 
the published structural catalogue. 

The simulations were also used to provide a reasonable faint 
limit in surface brightness that gives realistic results. A conserva- 
tive approach was used to define this limit, which we also set to be 
the same as that applied for the Gim2d fits. 



4.3 Small-size limit of GALFIT 

With the discovery of 16 new Ultra Compact Dwarfs in the Coma 
Cluster (Chiboucas et al. 2009) it is important to quantify how small 
an effective radius we can measure. To see how well Galfit can 
recover radii and magnitudes for small sources, 20000 further sim- 



ulations were carried out. The parameter space covered by this new 
set of simulations is presented in Table |4] 

We find from this new set of simulations that the recov- 
ered effective radii are unbiased, however for very small sources 
(i?e < 0.5 pixels) Galfit sometimes falls back to a hard-coded 
lower limit of i?e = 0.01 pix. This means that even for the per- 
fect conditions assumed in the simulations, the number of recov- 
ered sources that have around 0.5 pix will be fairly incomplete. 
It is very difficult for GALFIT to differentiate between a genuine 
point source and a small, yet extended, < 0.5 pixel source. 
These simulations assume perfect knowledge about the PSF, and 
model only the background contribution to the noise. Furthermore, 
the sources were injected on a relatively empty frame (visit 90). To 
see how well Galfit performs on real point sources, we inspect 
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Jvlagiiitude 


ocrsic index 


a 
P 






0.5 < n< 1.5 


-1.70 


0.036 




1.5 < n< 2.5 


-1.30 


0.020 




2.5 < n< 3.5 


-1.65 


0.035 




3.5 < n< 8.0 


-2.62 


0.079 


log, -Tie ^^vjrr / lV10Qt3i ) 


iSersic index 


P 


OL 




0.5 < n < 1.5 


-3.90 


0.12 




1.5 < n< 2.5 


-3.57 


0.11 




2.5 < n< 3.5 


-4.42 


0.15 




3.5 < n < 8.0 


-3.58 


0.11 


logn(GF/Model) 


Sersic index 




a 




0.5 < n< 1.5 


-4.09 


0.13 




1.5 < n< 2.5 


-3.17 


0.092 




2.5 < n < 3.5 


-3.01 


0.084 




3.5 < n< 8.0 


-2.77 


0.075 



700 



Table 3. Table of coefficients required to use Equation [5] to estimate the 
statistical errors on the total magnitude, Re and n, for various ranges of 
output n, for the Galfit fits. 




v90 
vl9 



1 2 3 4 5 

Effective radius [pixels] 

Figure 8. Histogram of Galfit Measurements of the Re for the sources 
found in visits 19 and 90. The filled histogram represents the sources in 
visit 90, i.e. a field with few Coma galaxies. The open histogram is the 
histogram of visit 19, with many Coma galaxies and globular clusters. 



Parameter 


Range 


Remark 


X 


400.0... 4000.0 




y 


400.0... 4000.0 




mag 


22.0... 26.5 




Re 


0.1. ..5.0 


Log 


n 


0.5... 6.0 


Log 


Me 


<25.5 




ell 


0...0.8 




pos 


0.0... 180.0 





Table 4. Parameters of simulated galaxies for the small radii simulations. 
Same comments as for Table|2]apply. i?eff and n are logarithmically spaced. 
As with Table[2]the limit in surface brightness is defined in terms of jie not 

(M>e 

the effective radii of sources in visit 19, the visit which covers the 
galaxy NGC4874. A large number of the sources in this visit are 
thought to be globular clusters (Peng et al. 2010), which should be 
unresolved at the distance of 100 Mpc. Fig. |8] shows a histogram 
of measured effective radii for this visit. It looks as if the distribu- 
tion of sources is shown by a powerlaw, plus an additional group of 
sources distributed around Re ^ 0.5 pix, with a standard deviation 
~ 0.3 pix. We conclude that Galfit output giving Re ^ 0.5 pix 
is likely to come from point sources. From the shape of the blue 
histogram in Fig. |8]we infer that Re > 1 pixel provides a robust 
lower limit for GALFIT i?e's. 



5 GOODNESS OF FIT INDICES 

Although the Sersic model is a well-known and tested fitting func- 
tion, real galaxies are more complicated. They present, among 
many other features, stellar bulges, star forming regions, AGNs, 
spiral arms, extended haloes, and central star clusters. 

In this section we define two complementary diagnostic in- 
dices, each designed to address in a different way the question of 
whether the Sersic model is an adequate fit, given the available data. 



or whether a more complicated function or extra components are 
required. These diagnostics are calculated for all fits. 

Following Blakeslee et al. (2006) we define the Residual Flux 
Fraction (RFF) as: 

- FLUXJSO ' 

where the summation is over all pixels within the I SCARE A of that 
particular object, |Resi | is the absolute value of residual image ob- 
tained by subtracting the best-fit model from the real galaxy image, 
and FLUX_I SO is the total flux within the I SOAREA, which can be 
taken from SEXTRACTOR. 

This index measures that part of the sum of the modulus of 
the pixels in the residual image which cannot be explained by the 
experimental error. The image variance is obtained from the usual 
CCD equation as: 

2 _ 2 I *^ 

C^image — Bkg ~^ ~ ' 

where S is the value of the model for that pixel, and g is the effec- 
tive gain, which can be found from the SEXTRACTOR parameters 
for the particular image. 

RFF quantifies the residual signal which cannot be explained 
by arguing that the fitting codes found a suboptimal minimum. It is 
best understood as a hypothesis testing procedure. 

If the real galaxy had a pure Sersic profile, both Gim2d and 
Galfit could find a model providing an exact fit to the galaxy. 
However, even in this optimal case, the errors associated with the 
readout noise and photon shot noise imply that the residual image 
will not be blank. In the case of independent errors, the properties 
of the residual image would be very similar to those of gaussian 
white noise, with a spatially varying a. The expectation value of 
the sum of the absolute values of these residuals is ^ 0.8 x (Jimage. 
Therefore, the expectation value of the numerator of the RFF is 0.0, 
should galaxies be pure Sersic models. Since the denominator is a 
normalization factor, the expected value of the RFF itself is 0.0 for 
such a model. Any positive or negative deviation from a pure Sersic 
model will increase the RFF. From our own visual experience with 
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this index, we find that fits with a RFF larger than 0.11 indicate 
that a more complex fit is required. This number was agreed after 
independent experiments made by CH and RG. 

The RFF diagnostic does not work well for objects with large 
I SCARE AS, and low (/i)e. In these cases, as both the galaxy and 
the model decay towards zero at large radii, the outer areas will 
dominate the RFF calculation, and even though the fit may have 
complicated and highly non-gaussian residuals at small radii, RFF 
will still be small. 

Hence we require a second diagnostic, which should be partic- 
ularly sensitive to the central residuals. This is calculated from the 
same set of pixels within the isophotal area of the target galaxies as 
RFF 

This complementary parameter is the Excess Variance Index 
(EVI). It is defined as: 



EVI : 



^ ^image ^ 



(9) 



where cfr is the root mean square of the residuals within the 
I SCARE A, and < crfmage > is the mean value of the square of 
the noise within the same area. EVI gives a measure of the gran- 
ularity excess of the actual residuals with respect to the expected 
granularity. 

For pure Sersic models, the numerator in the EVI expression 
above would simply be a sum of squares of normally distributed 
random variables with zero mean and pixel-dependent variance. 
The expectation value of the numerator would be, in this case, the 
value of the denominator. Hence, the parenthesized expression in 
the EVI definition has an expectation value of 0.0, and strong devi- 
ations from 0.0 indicate that the underlying real galaxy differs sig- 
nificantly from a Sersic model. Given the EVI definition, such de- 
viations would most likely be caused by the points with the largest 
intrinsic standard deviations (i.e., the inner points). This makes this 
index particularly sensitive to the structure of the residuals at small 
radii. The 1/3 factor in the definition of the EVI was later added 
so that residual images with an EVI larger than approximately 1.0 
present complicated substructure in the inner parts of the fits, whilst 
objects with an EVI smaller than approximately 1.0 show accept- 
able fits. In the end, after several iterations of this by eye calibration 
carried out by CH and RG, it was concluded that residuals with an 
EVI larger than 0.95 indicate complex residuals. 



5.1 Comparison of the diagnostic indices 

Since the RFF and EVI indices were designed to see whether a 
single-Sersic fit is sufficient, it is reasonable to ask how well the 
Gim2d and Galfit agree on this issue. As a first check, we plot 
in Fig.[9]the RFF and EVI indices of Galfit against their Gim2d 
counterparts. There is a strong correlation between the indices for 
the two codes. However, there is also a small offset present, in the 
sense that the Gim2d indices are on average slightly higher than 
the Galfit indices. We attribute this to a different treatment of the 
noise in the calculation of the two indices (the Galfit indices were 
calculated using the sky noise from the IVM maps, whereas for the 
Gim2d indices we used the noise as estimated by Gim2d itself.) 
For high values of RFF and EVI the correlation breaks down. Espe- 
cially for EVI, there are galaxies where the Gim2d index is normal, 
but the Galfit index is larger. In part this is the result of fitting low 
S/N sources, where Galfit has a preference for fitting a low sur- 
face brightness model to fluctuations in the background instead of 
fitting the source itself. 




0.1 0.2 0.3 

RFF(GALFIT) 




1 2 3 

EVI(GALFIT) 



Figure 9. Top panel: Residual Flux Fraction from Galfit plotted against 
that from Gim2d; Lower panel: Excess Variance Index from Galfit plot- 
ted against that from Gim2d. 



In Fig. [T0| we show histograms of the Residual Flux Frac- 
tion for different magnitude bins. If all sources were perfect Sersic 
galaxies, these should take the form of a Gaussian, where the width 
is dependent on the signal-to-noise. We find many faint sources 
with RFF < 0. We interpret this as overfitting: RFF < means 
that there is less noise in the image than expected, so, the code 
has modelled the noise away. For larger galaxies this problem is of 
course less severe, as the code has to model more pixels with higher 
signal-to-noise, so it has less freedom. 

Fig. pT] shows the calculated Excess Variance Indices, for the 
same magnitude bins. Although the faint sources centre around 
zero, the brightest magnitude bins show a large fraction of high 
EVIs. It is only in the brightest galaxies that we can see the devia- 
tion from the Sersic profile (top panels), because of the high signal- 
to-noise and spatial resolution. For the faint sources, the data are 
not good enough to detect any granularity. 

These indicators will allow separation into subsamples of 
galaxies which are well represented by the Sersic function, for stud- 
ies of the scaling relations or for separation by morphological class, 
or of galaxies for which additional components, truncation or ex- 
tension of the light distributions show that they must be fitted with 
more complex functions. 



6 THE CATALOGUE 

This section introduces the main result of the present paper: the 
catalogue of structural parameters for the Coma- ACS survey based 
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Figure 10. The Residual Flux Fraction calculated for our magnitude, size, 
and surface brightness clipped sample on the residual images of both codes. 
The sample is divided in three magnitude bins: 17.5-19.0 in the upper panel, 
19.0-22.0 in the middle panel and 22.0-25.0 in the lower panel. 
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Galfit Excess Variance Index G2D Excess Variance Index 

Figure 11. The Excess Variance Index, see caption of Fig. |10| for more 
explanation. 



on single Sersic fits. We provide a description of the catalogue and 
its limits in magnitude, radius, and effective surface brightness. 

The parameters presented in the final catalogue are summa- 
rized in Table [5] which presents the table headings of the released 
source lists together with a brief explanation of their meaning. In 
addition to best-fit values for both Gim2d and Galfit fitting, the 
catalogue includes coordinates and S EXTRACTOR photometry, and 
goodness of fit RFF and EVI for both codes. 

Recalling that all 74992 galaxies from the SEXTRACTOR pho- 
tometric catalogue from Paper II were fitted with Gim2d and Gal- 
fit, the full structural catalogue presented here contains all galaxies 
from the complete photometric catalogue which fulfill the selection 
criteria described in Sect. 16. II This amounts to 8814 sources. The 
full version of the structural parameter catalogue can be found in 
the electronic version of this paper, in the ASTRO-WISE website 
and the HST MAST archive. 
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COMA_ID 

RA ( J2 ) 
DEC (J2000) 
MAG_AUTO_CORR. 



F814W (GF) 

(7 (GF) 

Re (GF) 

CTR^ (GF) 

n (GF) 
(Jn (GF) 

Me (GF). 

(M>e (GF) 

Ellip. (GF) 

O'Ellip (GF) 

(GF) 

(70 (GF) 
RFF (GF) 

EVI (GF) 
F814W (G2D) 

0'F814W (G2D) 

Re (G2D) 

CTRe (G2D) 

n (G2D) 
CTn (G2D) 

Me (G2D) 

(M>e (G2D) 

Ellip. (G2D) 

O'Ellip (G2D) 
Pos Ang. (G2D) 

(TO (G2D) 
RFF (G2D) 

EVI (G2D) 



Name of source, as it appears in Paper 
II. 

Right ascension of source. 
Declination of source. 
Automatic F814W magnitude from 
SEXTRACTOR catalogue, corrected ac- 
cording to the prescription discussed in 
Paper II. 

F814W apparent magnitude of Galfit 
model. No k- or extinction correction 
applied. 

Error on F814W, as calculated from the 
simulations described in the text. 
Effective radius, in pixels from Galfit 
model. 

Error on Re, as calculated from the sim- 
ulations described in the text. 
Sersic index from Galfit model. 
Error on n, as calculated from the simu- 
lations described in the text. 
Surface brightness at Re calculated 
from the Galfit model. 
Mean surface brightness within Re cal- 
culated from the Galfit model. 
Ellipticity of model (l — b/a). Errors as 
given by Galfit. 

Formal error on ellipticity from Gal- 
fit. 

Position Angle of source. Errors as 
given by Galfit. 
Formal error on from Galfit. 
Residual Flux Fraction diagnostic, de- 
fined in the text, derived from the resid- 
uals from the Galfit model. 
Excess Variance Index, derived from the 
residuals from the Galfit model. 
F814W apparent magnitude of Gim2d 
model. No k- or extinction correction 
applied. 

Error on F814W, as calculated from the 
simulations described in the text. 
Effective radius, in pixels from Gim2d 
model. 

Error on Re, as calculated from the sim- 
ulations described in the text. 
Sersic index from Gim2d model. 
Error on n, as calculated from the simu- 
lations described in the text. 
Surface brightness at Re in the Gim2d 
model. 

Mean surface brightness within Re cal- 
culated from the Gim2d model. 
Ellipticity of model (l — b/a). Errors as 
yielded by Gim2d. 

Formal error on ellipticity from Gim2d. 
Position Angle of source. Errors as 
given by Gim2d. 
Formal error on from Gim2d. 
Residual Flux Fraction diagnostic, de- 
fined in the text, derived from the resid- 
uals from the Gim2d model. 
Excess Variance Index, derived from the 
residuals from the Gim2d model. 



Table 5. Output parameters included in the structural parameter catalogue. 
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6.1 Surface brightness, magnitude and size limits 

Although the two codes were used to fit all the sources detected by 
S Extractor in the F814W ACS images, not all of the resulting 
output parameters are meaningful and we apply limits in radius, 
magnitude and surface brightness to the final table, explained be- 
low. 

The Monte Carlo simulations presented in subsections §3.4| 
and §4.2| indicate that the uncertainty and reliability of the out- 
put parameters depend critically on the S/N ratio of the fitted 
sources. This is in turn a function of magnitude and surface bright- 
ness. Based upon these simulations we find that the derived pa- 
rameters and their errors are reliable for all sources for which: 
FSIAW < 24.5, and (/i)e < 24.5. 

Most galaxies for which F814VK ^ 24.5 are likely to be 
background galaxies, and in any case are beyond the magnitude 
limit of current spectroscopic surveys, so cluster membership can- 
not be verified, and their use in structural studies would be limited. 
At the distance of the Coma cluster, the apparent magnitude limit 
corresponds to Mf814w = —10.5, which is well within the abso- 
lute magnitude distribution of globular clusters. 

We do find, however, that a number of galaxies with compara- 
tively high central surface brightness, large Re, and in some cases 
with measured redshifts, fall below the surface brightness limit of 
(/i)e = 24.5. For this reason, we also include in our catalogue 
galaxies with 26.0 > (/i)e > 24.5, but we caution that because 
our simulations largely did not cover the parameter space occupied 
by these galaxies, we have less confidence in the derived structural 
parameters and their errors. For these sources for which only the 
central regions are detected, both codes are naturally forced to ex- 
trapolate a substantial part of the total surface brightness profile. In 
these cases, the results critically depend on the different hypothesis 
(e.g. constant sky) with which the codes work, and thus the derived 
parameters are more uncertain. 

In the final catalogue we include only sources for which 
FSIAW < 24.5 as measured with both codes, provided that both 
codes had converged. If only one code converged then the magni- 
tude from that code is used. For inclusion, sources have to satisfy 
the surface brightness criterion for either of the codes, not both. 

Besides magnitude and surface brightness cuts, we impose an 
Re lower limit with the aim to eliminate point sources from the 
catalogue, since for these the parameters of the Sersic fit have no 
meaning. From the simulations described in Sect. |4.3| we have de- 
cided to reject sources for which either code measures Re < 1.1 
pixels, which is 2a above the mean of the distribution in Figure [8] 
A number of point sources remain in the catalogue, in particular 
some bright, saturated stars whose wings give a larger Re. Follow- 
ing visual inspection of our images we estimate that < 2% of the 
objects in our sample are unresolved, these are a mixture of un- 
rejected stars, and unresolved objects such as Coma cluster Ultra- 
Compact dwarfs. 



6.2 Magnitude-surface brightness relation 

Although we explicitly exclude any physical analysis of the struc- 
tural catalogue in the present paper, we give a flavour of the types 
of objects included in the catalogue by showing (Fig.[T2| a plot of 
(/i)e against magnitude, measured with Galfit, for the objects in 
the catalogue. In this plot, blue denotes confirmed cluster members 
(from the redshifts of Marzke et al. 2010), red circles confirmed 
non-members, and the black dots objects without measured red- 
shifts. 
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Figure 12. (/i)e plotted against F814W magnitude for the objects in the 
catalogue. Blue points are spectroscopically confirmed cluster members, 
red circles are spectroscopically confirmed not to be members. Black dots 
without a red circle are objects with no measured redshift. The diagonal 
line of black dots at upper right is the remaining, partly saturated, stellar 
contamination. 



There is a diagonal line of black (and a few red) points at the 
upper right, this is the residual stellar contamination, in the case of 
the bright points these stars are saturated in the ACS images, and 
hence have measured Re > 1.1 pixels. The cluster members form 
a sequence towards the lower left, with a positive correlation be- 
tween (/i)e and luminosity. This sequence has considerable scatter 
as the galaxies have not been selected by morphological type nor n. 
Confirmed non-members (almost all background) have in general 
higher surface brightness for their apparent magnitude; however, 
there are also cluster members with brighter (/i)e. These are com- 
pact ellipticals similar to those discussed by Price et al. (2009). 



6.3 External Comparison 

Gutierrez et al. (2004) and Aguerri et al.(2005) study the struc- 
tural parameters of dwarf galaxies in the Coma cluster using 
ground-based data obtained with the Isaac Newton Telescope 
at La Palma, in the Sloan-r filter. Some of the decompositions 
in these papers fit multiple components, here we only focus on 
the galaxies which have bulge-to- total ratios B/T ^ 0.8 in the 
Gutierrez et al. sample and the dEs (excluding the dSOs) from 
the Aguerri et al. sample. 17 galaxies from Gutierrez et al. and 
7 galaxies from Aguerri et al. match within 2.0 arcseconds with 
galaxies in our structural catalogue. In Figure [13] we present a 
comparison between the effective radii Re and /Xe, the surface 
brightness at Re, for these galaxies. The correlation is good, 
with a few outliers, mostly galaxies with complex structures in 
their centres which are not well resolved in the ground-based 
data. 
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Figure 13. Comparison of our derived structural parameters with those of Aguerri et al. (2005) and Gutierrez et al. (2004). Left panel: Re in arcseconds; Right 
panel /ie- Filled symbols represent our Galfit values, and crosses our Gim2d values, as described in the legend at upper left in each panel. 



7 COMPARISON BETWEEN GIM2D AND GALFIT 

Both the Gim2d and Galfit simulations have shown that these 
two codes are capable of recovering parameters of simulated galax- 
ies in an almost unbiased way down to very low S/N values. How- 
ever, galaxies seldom consist of a single Sersic component. This 
makes a comparison between Galfit and Gim2d for real sources 
interesting. Even though both codes try to minimise a Figure-of- 
Merit (a value), there are significant differences in both the 
way this minimisation is achieved (Markov Chain Monte Carlo vs. 
Levenberg-Marquardt minimisation) and how the best fit parame- 
ters are determined (median value of a number of realizations in 
the case of Gim2d vs. final value after a number of iterations in the 
case of Galfit ). These differences could, in principle, have an im- 
pact on the performance of both codes even in the case of objects 
satisfying the magnitude, (/i)e and Re cuts described in §[6] 

7.1 Comparison of the fitted values 

The top two quadrants of Fig. [14] show the difference between 
the measured magnitudes, effective radii and Sersic indices plot- 
ted against the measured F814W magnitude. Quadrant (a) presents 
the results for objects with a measured Sersic index lower than 
2.5, whilst quadrant (b) presents the results for objects of higher 
Sersic indices. Quadrants (c) and (d) of this figure present these 
differences plotted against the F814W (/i)e for the same division 
in Sersic index. The objects included in these plots are a subset of 
the full catalogue described in § [6] The additional constraints for 
inclusion in these plots are: 

(i) (/i)e ^ 24.5 

(ii) RFF ^0.11 and EVI ^ 0.95 

(iii) 0.5 < n < 8.0 

In each case the galaxies had to satisfy the constraint for both 
codes. These additional cuts guarantee that the objects included in 
these plots are well measured and well represented by Sersic pro- 
files. 

These plots show that, whenever the two codes yield a Sersic 
index lower than 2.5, the observed agreement between the output 



parameters is quite remarkable even for the Sersic index itself. In 
this case, the expected magnitude scatter between the two codes is 
around 0.25 mag at the faintest levels, with a similar agreement in 
the Re measurements. 

Quadrants (b) and (d) of Fig. [14] show that in the cases in 
which the two codes measure a Sersic index higher than 2.5, the 
disagreement between the two codes is larger. It is therefore more 
difficult to get consistent measurements for objects with extended 
wings containing a large fraction of the object's flux. This disagree- 
ment between the two codes manifests itself not only as a larger 
scatter, but also as a systematic trend of the effective radii and 
Sersic index residuals, when plotted against the measured surface 
brightnesses. 

Whenever one code finds a brighter magnitude than the other, 
it also finds Re to be larger. This is in agreement with Fig. [3] This is 
probably caused by differences in the sky values adopted by the two 
codes. This is also the probable cause of the poorer agreement be- 
tween the two codes when fitting objects with high Sersic indices. 
For these objects, and in particular for sources that are found to be 
in the luminous haloes of the largest galaxies, it is very difficult 
to define the sky level with accuracy simply because these sources 
have a large fraction of their total flux in very extended wings and 
thus, any small discrepancy in the measured sky level translates into 
a large difference in the final luminosity and effective radius. 

Fig.[T4]also shows that the scatter in the output Sersic index is 
larger for higher values of this index. This is caused by the fact that 
models of high Sersic index start to converge towards a limiting 

profile and it is therefore very difficult to distinguish between 
two such models. 

Another discrepancy between the two codes comes from small 
sources. Often these sources are fit by Galfit with a very high 
Sersic index, and a half-light radius of a few pixels. These fits have 
a substantial portion of their light in wings which are below the 
detection threshold. On the other hand, Gim2d fits these sources 
with lower Sersic indices, and the Gim2d fit is probably a more 
realistic representation of the detected image. 

One advantage of Galfit is that it can model a number of 
galaxies simultaneously, leading to more consistent fits in cases 
where the target galaxy has near neighbours. Gim2d sometimes 
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Figure 14. Quadrant (a): Magnitude, effective radius and Sersic index residuals as a function of Galfit magnitude (left panels) and Gim2d magnitude (right 
panels). This quadrant shows results for galaxies with n < 2.5. Quadrant (b): same for objects with n > 2.5. Quadrant (c): Magnitude, effective radius and 
Sersic index residuals as a function of Galfit {/i}e (left panels) and Gim2d {/i}e (right panels). This quadrant shows results for galaxies with n < 2.5. 
Quadrant (d): same for objects with n > 2.5. Details of the vertical error bars and colour coding, are as in Fig.|4] 



has problems with these sources as the only thing Gim2d can do 
with neighbouring sources is to mask them, and that fraction of the 
flux from these neighbouring sources which falls outside of their 
masks will compromise the fits. 



8 SUMMARY, CONCLUSIONS AND FUTURE WORK 

We have fitted single Sersic models to the surface brightness dis- 
tributions of 74992 galaxies included in the photometric catalogue 
presented in Paper II, using the two most widely used two dimen- 
sional galaxy fitting codes, Gim2d and Galfit. Both codes create 



a PSF convolved trial function which is then compared to the in- 
put science data. The PSFs used were created mimicking the most 
important reduction steps that the science images had undergone. 

Independent simulations show that both codes can achieve 
similar accuracies for most cases. However, Gim2d requires a 
much higher degree of intervention compared with GALFIT to pro- 
duce the results. This is due to the greater flexibility of Galfit, 
which makes it easier for any given input image to satisfy the fit- 
ting hypothesis under which GALFIT works. Gim2d needs more 
manual work, including the use of enlarged masks with sizes cal- 
culated according to our detailed S extractor simulations, and the 
masking of nearby sources, to produce reliable results. 
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These simulations, together with an additional constraint 
aimed at rejecting probable point sources have been used to select 
a subsample of galaxies from the initial Sextractor catalogue. The 
criteria used are: 

• {/j)e < 26.0, where (/x)e is the mean surface brightness en- 
closed within Re. 

• F8UW < 24.5 

• Re > l.l pixels 
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We find that 8814 objects satisfy these criteria, cross-matching 
with the redshift catalogue of Marzke et al (2010), indicates that 
424 of these have redshifts, of which 163 have redshifts which place 
them as probable cluster members. 

We have also introduced two different goodness of fit diag- 
nostics. The Residual Elux Fraction is an indicator of the amount 
of light contained in the residuals which cannot be accounted for 
by the per-pixel photometric errors associated with the science im- 
age. The Excess Variance Index measures the importance of the in- 
ternal structure present in luminous galaxies that the Sersic model 
is unable to reproduce. These diagnostics indicate whether the sur- 
face brightness distributions of the fitted galaxies are well described 
by the Sersic profiles, or alternatively whether it is valuable to 
use more complex functional forms (maybe adding a nuclear point 
source, or a bulge+disk decomposition) than the simple Sersic pro- 
file used in this work. 

Although the Sersic fits presented here provide a good over- 
all description of the majority of galaxies in our sample, the spatial 
resolution and stability of the PSE provided by HST and ACS al- 
low the surface brightness to be measured accurately over a wide 
dynamic range. A range of more complex functions have been fit 
to surface brightness profiles, such as the "dual power law" (Fer- 
rarese et al. 1994), the "Nuker Law" (Lauer et al. 1995), and the 
core-Sersic profile (Graham et al. 2003; Trujillo et al. 2004). In 
the context of elliptical and spheroidal galaxies, a number of au- 
thors (e.g. Graham and Guzman 2003; Graham 2004; Ferrarese et 
al. 2006; Cote et al. 2007; Merritt & Milosavljevic 2005; Kormendy 
et al. 2009), have identified and measured either "extra light" or 
"missing light" with respect to the Sersic formalism in their pro- 
files. In future papers in this series we will investigate in detail 
these deviations from the Sersic profile, and the correlation with 
kinematic properties (Emsellem et al. 2007; Cappellari et al. 2007; 
Krajnovic et al. 2008, van Zee et al. 2004, Toloba et al. 2009). 

An alternative formalism which we will test for surface bright- 
ness profiles of ellipticals derives from the modified isothermal 
models of King (1962, 1966) and Wilson (1975), which are very 
successful in reproducing the surface brightness distributions of 
star clusters, and of some, but not all, elliptical galaxies. 

In the current paper we analyse only the radial surface bright- 
ness profile, but there is also information in the azimuthal distribu- 
tion, and we will undertake decompositions into disk, bulge, and 
nuclear components where appropriate (Andredakis et al. 1995; 
Balcells et al. 2007a,b), which will permit the construction of more 
meaningful structural parameter diagrams for the different physical 
components. 

We will also investigate asymmetric deviations from the ba- 
sic Sersic formulation of the light distribution, using techniques for 
parameterising galaxy surface brightness such as CAS (Conselice 
2003) and Gini/M20 (Lotz et al. 2004), and for detecting and quan- 
tifying bars (Marinova et al. 2010). 
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APPENDIX A: COMPARISON TO GIM2D AND GALFIT 
RESULTS BY THE GEMS COLLABORATION. 

Wanting to build upon the experience obtained by other teams 
in using GALFIT and GIM2D on HST/ACS data, we decided 
that the best starting point is the extensive GALFIT/GIM2D 



comparison by the GEMS collaboration (Haussler et al. 2007) 
on mock HST/ACS observations. The GEMS collaboration 
provides public access to the analysis and actual data on the web 
(http : //www.mpia . de/GEMS/f itting_paper . html). 
We used their BulgeOOOl and DiskOOOl simulated observa- 
tions. These contain artificial galaxies with pure Sersic profiles 
with n = 4 or n = 1 respectively, with appropriate Poisson 
noise, on top of a noise background created by putting together 
patches of real sky observations from different images. These 
simulated images are created for the F850LP filter using a PSF 
with a FWHM typical of this filter. The artificial galaxies cover a 
significant fraction of the galaxy parameter space expected for the 
Coma HST/ACS data, in terms of S/N ratio, total magnitude and 
effective radius. Therefore our analysis of these data is useful even 
though the Coma Legacy Survey does not use the F850LP filter. 

Following the work by Haussler et al. (2007) we fit a pure 
Sersic model with a flat background. All the parameters, including 
the Sersic index, were allowed to vary freely, although the value 
of the residual sky was later fixed, as explained below. The point 
spread function used was identical to the one used to create the 
mock galaxies. 

For Galfit we converged on a very similar fitting setup and 
best fitting results as Haussler et al.. For Gim2d we arrived at a 
different setup and masking treatment which led to improvements 
in the best fitting results. 

First, Gim2d was run with its default decision algorithm on 
the size of the cutout to be fitted. It was also allowed to fit and refine 
the sky value by itself (DOBKG=YES) and to automatically estimate 
initial values for the fitted parameters (INITPARAM=YES). The 
Sersic index was left as a free floating parameter. The S Extrac- 
tor parameter configuration ensures detection of sources with a 
global signal-to-noise ratio of 20.25. The SEXTRACTOR parameter 
BACK SIZE was set to 512. The minimum stamp size was set to 
31px. This experimental configuration leads to the "recommended" 
setup for Gim2d. The noise model used for the Gim2d fits was 
based on the traditional CCD equation using the background a and 
the image effective gain which is ultimately regulated by the expo- 
sure time and reduction process. 

For the DiskOOOl models Gim2d produced a fit for 470 sim- 
ulated profiles. Figure [aT] shows the results obtained by Gim2d in 
this case. The upper panel presents the magnitude residuals, the 
middle panel presents the ratio in effective radii, and the lower 
panel presents the output Sersic index as a function of the average 
surface brightness within a effective radius of the simulation. Our 
Gim2d results are shown by red crosses and the GEMS team re- 
sults with Gim2d by the blue circles. The general conclusion from 
this is that there are systematic trends at the faintest surface bright- 
ness levels, fainter than (/i)e ^ 22.5. The scatter between fit and 
simulation is smaller for our Gim2d run compared with the results 
of the GEMS team. Likely,this is caused by the fact that we masked 
overlapping sources with a value of -2 in the mask image, a situ- 
ation that forces Gim2d to ignore those problematic pixels. The 
GEMS group did not do this. 

For the BulgeOOOl models Gim2d yielded a successful fit for 
558 simulations. Figure [A2] shows the results. The general conclu- 
sion is that our results are an improvement on the GEMS results 
down to a surface brightness of around z ^ 22.5, and have a 
smaller scatter. However there are again systematic trends which 
increase with fainter surface brightness. 

We investigated if the systematic deviations shown by the 
Gim2d solutions at surface brightness levels fainter than z — 23.0 
could be due to an inaccurate sky. To this end a number of low 
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Figure Al. Results for GIM2D on pure exponential models. From top to 
bottom: magnitude residual (GIM2D minus Model), ratio of effective radius 
(Gim2d / Model) and Sersic index as a function of the true average surface 
brightness within an effective radius. Red crosses are our GIMD2D runs, 
blue crosses are the results for the same simulated galaxies from the GEMS 
team. 



Figure A2. Results for GIM2D on pure de Vaucouleurs models. From top 
to bottom: magnitude residual (GIM2D minus Model), ratio of effective 
radius (Gim2d / Model) and Sersic index as a function of the true average 
surface brightness within an effective radius. Red crosses are our GIMD2D 
runs, blue crosses are the results for the same simulated galaxies from the 
GEMS team. 



surface brightness sources were refit. This time, we used a sky de- 
termination derived from a customized configuration file for S EX- 
TRACTOR whose main difference with the S EXTRACTOR config- 
uration file used in the previous fitting round is that it uses a 2048 
mesh for background subtraction and determination. Also, we did 
not allow Gim2d to fit the sky. Furthermore the cutout size for 
Gim2d was increased to 25 times the default size. 

Figure |A3| shows the result of this test for the mock galaxies 
present in image DiskOOOl. Comparison with Figure |A1| shows 
that, using this prescription to deal with the sky level, Gim2d no 
longer shows systematic offsets in magnitudes and effective radius 
even at the faintest surface brightness levels. 

Figure |A4| shows the result of this fitting round for the mod- 
els included in the BulgeOOOl image. Comparison to Figure |A2] 
clearly shows that setting the sky to be equal to the sky value given 
by S Extractor makes the systematic offsets in magnitudes and 
effective radius to vanish even at the faintest surface brightness lev- 
els. The statistical errors (the scatter) are rather large however. 

From the fits to the synthetic models generated by the GEMS 
team comprising either pure exponentials or pure de Vaucouleurs 
models, we conclude that: 

• The errors in the recovery of parameters appear to be domi- 
nated by errors in the sky determination. 

• Our G2D fits show significantly reduced scatter with respect 
to the G2D fits published by GEMS team. The difference can prob- 
ably be traced to the manual masking of intervening sources in the 
object's mask done by our team. 

• Increasing the cutout size used by Gim2d reduces or removes 
the systematic offsets in magnitude and effective radius found when 
the default cutout size is used. 

• Overall, we conclude that we can use both the GALFIT and 
GIM2D codes with a similar or better accuracy to that found by the 
GEMS collaboration. 
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Figure A3. Results for GIM2D on pure exponential models of low surface 
brightness, using our own sky determination (see text). From top to bot- 
tom: magnitude residual (GIM2D minus Model), ratio of effective radius 
(Gim2d / Model) and Sersic index as a function of the true average surface 
brightness within an effective radius. Red crosses are our GIMD2D runs, 
blue crosses are the results for the same simulated galaxies from the GEMS 
team. 



APPENDIX B: CREATING OBJECT MASKS FOR GIM2D 

We provide a detailed recipe for producing object masks for Gim2d 
that solves the missing flux problem that appears when standard 
S Extractor segmentation masks are used in Gim2d. 

(i) The aperture for the target should ideally be equal to the 
Kron aperture of an object with the same FLUX_RADIUS and n 
that would be measured if the observations had infinite S/N ratio. 
The elongation and position angle of this aperture would ideally 
be equal to those of the real object. The pixels within this aper- 
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Figure A4. Results for GIM2D on pure de Vaucouleurs models of low sur- 
face brightness, using our own sky determination (see text). From top to 
bottom: magnitude residual (GIM2D minus Model), ratio of effective ra- 
dius (Gim2d / Model) and Sersic index as a function of the true average 
surface brightness within an effective radius. Red crosses are our GIMD2D 
runs, blue crosses are the results for the same simulated galaxies from the 
GEMS team. 



ground sources are flagged as 0, so that Gim2d regards these pixels 
as sky. 

In a few cases, for instance, with objects with close neigh- 
bours or near to the CCD edges, S Extractor did not provide 
measurements of MAG_AUTO and FLUX_RADIUS. In these cases 
MUOBS is estimated by spreading the light encircled in the isopho- 
tal area of the object (ISOAREA_IMAGE) uniformly in an aperture 
with a radius equal to | x ^ isQARea_im age ^ j^^^ ^^^^^ 
be somewhat brighter than a measured MUOBS, although the final 
apertures obtained were similar in size. 

The size of the poststage stamp over which Gim2d has to 
work is an integral part of the suggested solution. The quantity 
FRAD^OD amounts, in the vast majority of the cases, to 50% 
of FLUX_RADIUS, as calculated by SEXTRACTOR. For the ob- 
jects where SEXTRACTOR could not calculate FLUX_RADIUS, 
FRAD^OD is half of the radius of the circle with an area equal 
to ISOAREA_IMAGE, which always does exist. 

The number YFX (MUOBS) is given by the empirical rela- 
tion between the 50% FLUX_RADIUS and the input effective ra- 
dius that is appropriate for sources with a high n. For reasonable 
Sersic indices YFX (MUOBS) is given by the fifth-order polyno- 
mial in the mean observed surface brightness MUOBS, presented in 
paper II. The final size of the post stamps imagelets extracted is 
7.0 X FRAD_MOD x YFX(MUOBS) on a side. 



ture should be flagged as 1 in the mask file so that Gim2d recog- 
nizes them as belonging to its target object. This aperture is cal- 
culated from the SEXTRACTOR measurements. FLUX_RADIUS, 
THETA.IMAGE and ELONGATION were taken directly from the 
SEXTRACTOR catalogue. MUOBS is a SEXTRACTOR based esti- 
mate of effective surface brightness of that particular detection, ob- 
tained as if MAG_AUTO were the total magnitude and the 50% of 
FLUX_RADIUS were the real effective radius. It is thus an esti- 
mate of the effective surface brightness. The 5th degree polyno- 
mials used in the SEXTRACTOR simulations presented in Paper II 
were then used to estimate empirically the effective radius of the 
source. This effective radius estimate is then converted into a Kron 
aperture using the relationship between Re and Ri (the S/N = oo 
Kron aperture in units of the effective radius) and n. n = 2.5 is 
assumed where n is undefined. Although the aperture defined in 
this way leaves out a finite fraction of the total flux from the target 
models outside the mask, this magnitude offset is 0. 1 mag in the 
worst case. For this reason, the initial flux estimate will be around 
0.1 mag fainter than the true integrated flux. As above, a very simi- 
lar statement can be made about the initial effective radii estimates. 
Using this customized mask, Gim2d has a realistic opportunity to 
measure 100% of the target object total flux and size, whereas with 
the default SEXTRACTOR segmentation image, this is not true. 

(ii) The apertures for the background sources are elliptical aper- 
tures, with an area given by the ISOAREA.IMAGE SEXTRACTOR 
measurement. The elongation and position angles of these ellipses 
are also taken from the SEXTRACTOR measurements. These pixels 
are flagged as -2 so that Gim2d ignores these pixels in all cal- 
culations. Gim2d is therefore blind to the innermost regions of the 
background sources. Thus only the outer regions of the background 
sources contribute to the sky level affecting the galaxy under study. 
In the case of an overlap between target and background object pix- 
els, the pixels are assigned to the background object as including 
them in the target would potentially compromise the fit. 

(iii) Finally, pixels not belonging to the target or to the back- 



